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ABSTRACT 


The ionospheric response (i.e., the sensitivity and delay of the 
electron content) to the 270day fluctuations of solar radiation has been 
theoretically analyzed. In the context of the present work the sensi- 
tivity is defined as the magnitude of the electron content fluctuation 
for a given change in solar flux, and the delay is defined as the time- 
shift of the response of the electron content to fluctuations in solar 
flux. Both the model in which the action of neutral winds on the F-2 
layer is disregarded and that in which the wind effect is included as 
part of a positive feedback mechanism, are studied. It is shown that 
the neutral winds play a dominant role in the mentioned ionospheric 
response. The computed delays decrease with increasing solar activity 
in both models. Nevertheless, when the action of neutral winds is con- 
sidered the delays are larger. The experimental verification of the 
delays computed theoretically is only possible when two conditions are 
met: First, that the solar flux fluctuations exhibit sufficient ampli- 

tude, and secondly, that the first condition occur during geomagnetically 
quiet periods. For the level of solar activity exhibited during the 
year period when observations were made, a fairly good agreement was 
found, during the favorable period in which the two conditions mentioned 
were met, between the delay predicted by the model in which neutral winds 
are included and the experimental results. The computed sensitivities 
from this model have found a semiquantitative agreement with experimental 
results due to the presence of other effects not included in the model. 
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CHAPTER 1 


INTRODUCTION 

One of the main processes taking place in the ionosphere is the 
photoionization of the neutral atmosphere by solar extreme ultraviolet 
radiation (EUV) . Most of the early papers concerning the influence of 
solar radiation on the ionosphere dealt with variations of the critical 
frequency, that is, of the maximum electron concentration, and the vir- 
tual height of the layers obtained by vertical incidence sounding at 
numerous observatories over a period of many years. The values of the 
critical frequency and virtual height are derived from examination of 
ionograms, which though very rich in information, are of difficult 
interpretation and very expensive in their processing. 

It was realized very soon that virtual height may be a misleading 
parameter in that the time of travel of a radio pulse can be greatly 
influenced by changes in the distribution of electrons well below the 
level of reflection as well as by changes in height of the reflection 
point itself. Consequently, all efforts were focused on critical fre- 
quencies studies. The chief result of that epoch can be synthesized 
by the expression: 

f _F2 oc N F2 cc 1 + 0.02R 
0 m 

where R is the mean Z\irich sunspot number, fgF2 is the critical frequency 
of the F2 layer and N F2, the maximum electron concentration. N F2 and 
fgF2 values correspond to any midlatitude station, at noontime 
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(Ratcliffe and Weekes[l960]) . 


The idea that the above expression could represent the worldwide 

results has long since been discarded. In fact, at midlatitude stations 

noon values of N F2 are greater in winter than in summer (seasonal 
m 

anomaly). On the other hand, studies of worldwide f Q F2 data have shown 
that the anomaly is not exclusively seasonal but contains annual and 
semiannual components also ( Rishbeth [l968 ]) . 

With the launching of artificial satellites it was possible to 
record in ground based stations information about the columnar content 
of plasma in the ionosphere with excellent temporal resolution and very 
inexpensive equipment. In recent years the vast amount of observational 
electron content data, obtained mostly from measurements of the Faraday 
rotation angle of VHF telemetry transmissions from geostationary 
satellites, have proven to be a very valuable source of information in 
studying general aspects of ionospheric physics and transionospheric 
propagation, ( Garriott, da Rosa, and Ross [l970]), and in particular it 
is an excellent index of the ionospheric response to the solar activity. 

Bhonsle, da Rosa, and Garriott [1965] , using electron content infor- 
mation from 1958 to 1962 derived from the observation of low orbiting 
satellites, plotted the noontime values of electron content versus the 
smoothed sunspot number and observed that the data grouped themselves 
around three straight lines (see Fig. 1.1) of progressively smaller 
slopes. The maximum slope corresponded to the straight line connecting 
all equinoctial points, then followed the slope corresponding to the 
line connecting the winter points, and finally the smallest slope 
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Fig. 1.1. Solar cycle variations of I of midlati- 
tude ionosphere and the mean sunspot number, R. 
(After Bhonsle, da Rosa and Garriott[l965]> . 


corresponded to the straight line connecting all summer points. From 
the graph, they derived the empirical formula: 


I = [a + b(R - 40)] X 10 17 el.m 2 


(for R s 40) 


The parameters a and b are constants which are different for differ- 
ent seasons. The data used for these calculations were insufficient to 
permit a careful study of the seasonal variation of the above relation- 
ship and, in view of the relative motion between satellite and observer, 
these data did not correspond to a fixed ionospheric point. Consequently, 
all conclusions drawn from that study should be taken with some reser- 
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vations . 


da Rosa et al.[l972] in a preliminary investigation of the influ- 
ence of the solar activity on the ionospheric electron content used 
electron content, data from 1964 to 1970, obtained from the measurements 
of the Faraday rotation of VHF telemetry transmissions from ATS 
(Applications Technology Satellite) and SYNCOM geostationary satellites. 
In view of the lack of continuous monitoring of EUV radiation, and 
after an exhaustive cross-correlation analysis between four different 
solar activity parameters and several ionospheric quantities they 
adopted for the investigation of long-term effects the half-year mean 
of the solar radio flux at 2.8 GHz as the index of solar activity and 
the 31-day mean value of the daily average electron content as the in- 
dex of electron content. 

They concluded that at least over the range of the solar activity 
represented by the semi-cycle covered, the relationship between electron 
content and the long-term average of solar microwave flux can be ex- 
pressed by: 

IMEAN 31 = I Q + a(S> 183 

where cr is defined as the "sensitivity” of the ionosphere to long-term 
changes in solar activity. 

By assuming that long and short term sensitivities are identical 
it was possible to especialize the above formula to express the prog- 
nosticated values of the daily averages of the total electron content 
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corresponding to observed daily averages of solar radio flux, i.e.: 


IMEAN = + aS 

P 0 

Plots of both IMEANp and IMEAN, that is, of the prognosticated and 
of the observed values of mean daily electron content, showed that the 
above empirical relation seems to fit reasonably well the data for both 
low and high solar activity, although, an apparent time-shift between 

IMEAN and IMEAN was observed which suggests that the fluctuations in 

P 

the real ionosphere lag the corresponding fluctuations of solar flux. 

The present work deals with the ionospheric response to fast fluc- 
tuations in solar radiation such as those associated with a solar rota- 
tion. It is well known that the emission of solar radiation does not 
occur uniformly over the surface of the sun. Instead, it may be con- 
centrated on certain active areas that can persist for several months, 
particularly during the years of high solar activity. Consequently, as 
the sun rotates with a period of about 27 days, solar activity as seen 
from the earth will seem to vary with that period. The variation in 
the incident solar radiation is responsible for the appearance of 27- 
day variations in the earth’s ionosphere. In the absence of other 
effects, the total electron content during one of these 27-day periods 
will also vary accordingly and so do the rest of ionospheric parameters 
depending directly upon the photoionization process. 

This is basically the idea developed in Chapter 2 where, after 
showing that the daily averages of the variables involved in the conti- 
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nuity equation can be expressed in terms of the rotation frequency of 
the sun only, the continuity equation is solved in the absence of all 
kinds of drifts except those produced by ambipolar diffusion, to find 
the time-shift between the daily averages of the solar input and of the 
total electron content. It is found that the ionospheric delay de- 
creases when solar activity increases. The results found in this chap- 
ter do not correspond to the reality because they come from an over- 
simplified ionospheric model that has resulted from disregarding the 
action of neutral winds on the electron concentration profile. Never- 
theless they are relevant when compared with the results found in 
Chapter 4, where those drifts were considered, to point out the impor- 
tant role that the neutral winds have in the ionospheric response to 
the solar input. 

In Chapter 3, by using the general feature which describes the 
behavior of the neutral wind velocities, a feedback mechanism is pro- 
posed as part of a non-linear system which simulates the ionosphere. 

The input to that system is a series of step functions used to approxi- 
mate a sinusoidal wave of 27-days period representing the daily aver- 
ages of the solar flux. The response at the output of the system is 
the electron content. A detailed analysis of the equations involved 
in the system concludes this chapter. 

Chapter 4 has been divided into two parts: the first contains a 

description of the numerical method used for integrating the set of 
differential equations involved in the system discussed in Chapter 3, 
and the last presents the results obtained for periods of time corre- 
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sponding to the autumnal equinox of low, medium, and high solar activi- 
ties. 

The results found in Chapter 2 show that the sensitivity of the 
electron content, i.e., the relative variation of the electron content 
throughout a given period of time, is smaller than the corresponding 
relative variation of the solar flux, and that the delay of the response 
electron content to the solar flux decreases when solar activity in- 
creases. The conclusion concerning the sensitivity of the electron con- 
tent is modified later by the results obtained in Chapter 4, when the 
neutral winds effect on the ionospheric F2 region is included, and the 
conclusion concerning delay is emphasized. 

In Chapter 5, electron content data corresponding to eight differ- 
ent stations and during the period of one year, are statistically an- 
alyzed. A creiterion is established in terms of the total signal-to- 
noise ratio which allows for detecting those periods when the 27-days 
effect dominates, and the measured delays can be compared with the 
theoretical ones computed in Chapter 4. When this is done, the agree- 
ment is fairly good. 
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CHAPTER 2 


DERIVATION OF THE IONOSPHERIC TIME-RESPONSE TO FLUCTUATIONS IN SOLAR 
ACTIVITY FROM AN IDEALIZED SOLUTION OF THE CONTINUITY EQUATION 

In this chapter the equation of continuity is solved for the elec- 
tron content to determine the time delay between daily averages of elec- 
tron content and of the fluctuations of the solar radio flux caused by 
the 27-day changes in solar radiation. 

Let us consider the path between a ground based station placed at 
a midlatitude point in the northern hemisphere and a geostationary 
satellite, orbiting at 6.6 Earth radii (see Fig. 2.1). Assuming a 
linear loss coefficient over the whole path, the continuity equation 
for electrons is: 


3N _ 

at “ 


q 


$N - V.(NV) 


( 2 . 1 ) 


where N is the electron concentration 

q is the rate of ionization-production 
3 is the linear loss coefficient 


and V is the drift velocity due to ambipolar diffusion, neutral winds, 
and electric fields. 

(N, q, P, and 'V are functions of height and time.) 

Integrating along the path between the station (observer) and the 
satellite and calling 


I 



(Electron content) 
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sat . 



(Integrated production 
rate) 


(Weighted mean recombi- 
nation rate) 


Equation 2.1 becomes: 


sat . 


1 1 = Q ~ 0>I ~ j [V-(NV)]ds 
obs . 


( 2 . 2 ) 


In order to evaluate the relative importance of the integrated 
transport term (third term on the r.h.s. in Equation 2.2), it is useful 
to think of the electron content as divided into two parts: the iono- 

spheric content; and the protonospheric part. 

To a great extent, the value of the ionospheric electron content 
is determined by the content of the F2-region, ( da Rosa et al, op.cit. 
[1972]), (a layer consisting mainly of 0 + ions and electrons, the top 
of which can be placed at roughly 1000 km, where H + ions begin to be 
present in greater concentrations) ( Whitten and Poppof [1971]) . In this 
region the ionization is constrained to move along the geomagnetic 
field lines, except for electromagnetic drifts set up by transverse 
electric fields important only during storms ( Abur-Robb [l969] , Amayenc 
and Vasseur[l972] , Bramley and Young[1968], Stubbe and Chandra[l970]) . 
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Fig. 2.1. Sketch of a magnetic-meridian cross-section of the earth 
and its plasma environment, showing that for midlatitude points the 
lines of force of the earth's magnetic field are practically parallel 
throughout the ionosphere, which extends to roughly 1000 km. 

For midlatitudes, in the ionospheric range of heights, the lines 
of force of the geomagnetic field are practically parallel, as shown 
in Fig. 2.1, therefore, the divergence term can be approximated by: 


v -(™) = | e (N V 

where i is measured along the magnetic field line. 
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If is the value of s at the top of the ionosphere, the inte- 
grated transport term can be expressed: 


/ [v-On?)]ds -j 1 Ijdnr,) |f 61 


3s 

But = cos B ~ constant 

where 0 is the angle between the magnetic field line and the ray path, 
Therefore: 


c 

/ 


[V* (NV) ]ds 


= COS0 


J,. 

r 




)dj£ = cos0 


Z 1 

J 


d(NV^) 


^cosB 


where is the flux along the geomagnetic field line at the top of the 
ionosphere, i.e. at 1000 km height. 

In terms of the topside flux 4^, Equation 2.2 becomes; 

|£ = Q - <B>I - ^cosB (2.3) 


where all quantities now refer to the ionospheric F2-region. 

Whistler observations near L = 4 under quiet conditions have shown 

8 “2 —1 

a daytime upward flux of 3 X 10 el.m sec and a nighttime down- 

8 -2 “1 

ward flux $ of 1.5 X 10 el.m sec (Park[1970]) . Other observations 
i 

by incoherent scatter radar ( Evans [l971] , Carpenter and Bowhill [l97l]) 
are in substantial agreement with the whistler results. However, at 
midlatitudes, still under quiet conditions, whistler observations in- 
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dicate that the total interchange of flux between the ionosphere and 


the protonosphere in a 24 hour period is nearly zero ( Park, C.G. , 
personal communication). Therefore, if the equation of continuity as 
given by Equation 2.3 is especialized for daily averages, the last term 
of the r.h.s vanishes, leaving: 

|| = q - <pyr (2.4) 

On the other hand, when protonospheric content is considered, the 
transport term cannot be neglected in the continuity equation; in fact, 
it is the dominant cause of changes in protonospheric content since 
production and losses are practically nonexistent in that region. How- 
ever, if daily averages of protonospheric and ionospheric content 
changes are compared, the former can be disregarded based on the follow- 
ing reason: during daytime the protonospheric content amounts to a 

small fraction of the total electron content (it has been estimated 
at about 10% by da Rosa and Garriott [l969], Almeida [l972]). During 
the night, that fraction becomes important approaching some 50% of the 
total, but the ratio between daytime and nighttime total electron con- 
tent generally falls between 5 and 10. Therefore, on a daily average 
basis the protonospheric content represents a small fraction of the 
ionospheric content, and the meaning of Equation 2.4 can be extended to 
represent the general behavior of the quantities between the observer 
and the satellite. 
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In Equation 2.3 each of the quantities (which we will represent 
by the generic symbol x(t)) shows a diurnal variation that will be 
approximated by a sinusoidal law: 


x(t) = a Q + a^ cosa^t + a^ sinu^t 


(2.5) 


where U) E is the Earth's rotation frequency, i.e. uu^, = 2it rad. day \ 

In addition, these quantities show a day-to-day change of erratic 
nature, the cause of which is not yet completely understood. Superim- 
posed on such a random variation there is a systematic response to the 
fluctuations of solar ionizing energy. These fluctuations are, during 
periods of high solar activity, associated with discrete active areas 
in the sun and tend to behave in a sinusoidal fashion with a period of 
27 days (the synodical period of the sun’s rotation). Since we are 
investigating the ionospheric response to changes in solar EUV we will 
treat the erratic day-to-day changes as noise and consider the EUV 
fluctuations sinusoidal. Thus the coefficients of Equation 2.5 are 
themselves time-dependent: 


a.(t) = a._ + a. , cosUU t + a.„ sinuu t (i = 0,1,2) (2.6) 

l lO ll s i2 s ’ 


where UJ is the sun’s rotation frequency, i.e., u) = (l) / 27 . 

S S E 

The functions x(t) are then amplitude modulated sine functions, 
the modulation being sinusoidal. 
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Combining Equation 2,6 with 2.5 one obtains: 


x(t) = X Q + costi) g t + X 2 sina) g t + Xg cosuo^t + X^ sintt^t 

+ X cos(m + 0) ) t + X„ sin(uo + utf )t 
5 E s 6 E s 

+ X ? cos(u) E - UJ g )t + Xg sinCujg - uu g )t (2.7) 


where: X_. = a 


0 00 


X, = a 


01 


X„ = a 


02 


X„ = a 


10 


X, = a 


20 


X 5 = <a U - a 22>/2 
X 6 = <S 21 + a i2 ) / 2 

X 7 = (a ll + a 22 ) / 2 

X 8 * <a 21 - a i 2>/ 2 


The daily average of x(t) is, by definition: 


x(t) = 



( 2 . 8 ) 


Substituting the value of x(t), given by Equation (2.7), into the above 
integral and considering that: 


y 


T 

4- E 

t+ - 


T 

"E " _E 
2 


cos 


V 


dt 


-kJ 


t t e 
t+ - 


_ sinil) t dt 

E - V E 
2 


= 0 


(2.9) 
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cosuj t dt = 0.9977 cosuu t 
s s 


( 2 . 10 ) 


and 



sinii) t dt = 0.9977 sinuo t 
s s 


( 2 . 11 ) 



cos(U) ± uj )t dt = 0.017 cos(ti) ± UJ )t 
E s E s 


sin(UJ ± (J) )t dt = 0.017 sin(uu ± UJ )t 
E s E s 


( 2 . 12 ) 


(2.13) 


It can be seen that in view of the fact that the various coefficients 
are either near zero or near unity, the daily average of x(t) can be 
approximated by: 


x(t) = X n + X cosuu t + X„ sinU) t . (2.14) 

0 1 s 2 s 

Similarly, if we calculate first the time-derivative of both sides of 
Equation (2.7) and then its daily average, the result can be expressed 
by the equality: 

= 1-^(7) (2.15) 
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Finally, if we consider two functions x(t) and y(t), the diurnal 
variations of which are given by expressions of the form (2.6), the 
daily average of the product is: 


x(t)y(t) = 


kl 


T 

. E 
t+ ~ 

T ( a o +a i cosuo £ t +a 2 sinuUgt) (b^+b^ cosu) E t+b 2 sinuu E t)dt 
E 


where all coefficients on the r.h.s, of the equation present a time- 
dependence of the form given by Equation (2.6). After performing the 
product indicated under the integral sign, the r.h.s. of the equation 
is decomposed in a sum of integrals the values of which are given by 
the expressions labelled from Equation (2.9) to (2.13) from which the 
nearly-zero coefficients can be disregarded, yielding: 

x(t)y(t) = x(t) X y(t) + ^ (a.^ + a 2 b 2^ (2.16) 


Applying the results found to the equation of continuity as given by 
Equation (2.4), one obtains: 

ft 1 = Q - <P> X I - | [<0)^ + <0> 2 i 2 l (2.17) 

where the term between brackets depends directly on the amplitudes of 
the diurnal variations of (0) (weighted recombination rate) and of I 
(electron content); i.e.: 
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*, ■ (h 2 * t ) 1/2 

\ s >=(< s >? + < e > 2 ) 1/2 


(Amplitude of the diurnal variation of I) 
(Amplitude of the diurnal variation of (P)) 


Calling Tj the ratio between the last two terms on the r.h.s. of (2.17), 
i .e . : 


*n = 


<PVl + <P>2 1 2 

2<P> X I 



(2.18) 


where >is the phase difference between the diurnal maximum of (p) 

\P > 

and I, Equation (2.17) takes the form: 


i| ¥ T + <B>XI =|Q 


(2.19) 


where § = 1 + T] 

An approximate value of 7| was obtained using a simulation program 
of the ionosphere based on the coupled continuity and heat flow equations 
( Waldman [l971 ]) . The program which makes use of observed values of 
electron content as one of the necessary boundary conditions, can be 
operated both in a steady state and in a time-varying mode and its self- 
consistency is only limited by the fact that the neutral wind is entered 
as an input, not solved for in a closed form. The diurnal variations 
of (g) and I obtained from this program when especialized for a mid- 
latitude point, a solar radio flux of 175 units and zero wind velocity, 
are represented in Fig. (2.2). A sinusoidal curve was fitted through the 
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computed points (dash lines) and from it the ratios below were obtained: 


= 0.2947 

P 


= 0.85 
I 


and 


COS$^py= 0.81 


Substituting these values into (2.18), a value of 0.10 was found for T], 
On the other hand, if (2.19) is written symbolically in the 
following way: 


the first of the two terms between brackets is about three orders of 
magnitude smaller than the second, since: 


1 3 


uu 

s 


5 at - T 


TO ' 6 - 1 

10 sec 


and ( 0 ) ~ 10 ^ sec 


therefore Equation (2.19) becomes: 


<3> x I = i Q 


( 2 . 20 ) 
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LOCAL TIME (hours) 


J 

Fig. 2.2. Diurnal variations of ({3) and I, with their corresponding 
best-fit curves, obtained from a simulation program of the ionosphere 
based on the coupled continuity and heat flow equations (Waldman[l971 ]) 
for the case in which a solar radio flux of 175 units is used. 
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where (j3), I, and Q have the form given by (2.14); i.e.: 


<3> = (P >Q + <P> X cos^t + (p> 2 sinii^t 

I = I ft + I, COSUJ t + I sinu) t 
0 1 s 2 s 

Q = Q 0 + Qi coso) s t 

(Q has been conventionally taken with zero phase) . 

Substituting the expressions for (P) , I, and Q into (2.20), that 
equation becomes: 


GV ^)i cosu, w t + <3> 2 sinuL) s t )( I 0 +I i cost V +I 2 Simi) s t ) " f ( Q 0 +Q l COStJU s t ) = 0 

Since this - equation must apply at all values of time, the coefficients 
of cosu^t, sinu^t, and the time independent part must, separately * be 
equal to zero. Therefore the above equation is separated into three 
equations independent of time from which the following expression for 
the delay-time between I and Q is obtained: 


where 


6, = — arctan 
I UU 

s 



t 1 - 1 c ° sl " s ^>) slM 4e> (2-21) 

f e} + 3 sln V <E1 >- cosl "s 6 <e> 

(Ratio between the amplitude of the a.c. term 
and the mean, for production) 



(Ratio between the amplitude of the a.c. 
term and the mean, for recombination) 
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tan U) 6 
w 


< S >2 
<f>> = <e>! 


(Phase of the recombination rate with respect 
phase of the production term) 


and 


tanuu 6 t 
s I 



(Phase of the electron content with respect 
phase of the production term) 


The unknowns on the r.h.s. of Equation (2.21) are D^^and 6 
the determination of which in terms of the solar input can be carried 
out in the following way: 

By definition, the daily average of the weighted mean of the re- 
combination coefficient is: 




where m is the number of (0) values over which the mean is taken, 

0 is the linear loss coefficient (function of height and time), 
N is the electron concentration (function of height and time) , 
h Q is the lower boundary of the F2-region, 
and h is the independent variable height. 


In order to determine the r.h.s. of Equation (2.22) it is necessary 

to know the neutral composition of the atmosphere, since 0 is a linear 

combination of the neutral concentrations of N and 0 , and the electron 

Z z 

concentration profile. 
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For each individual gas in the atmosphere, the variation of its 

' / 

concentration with height is given by the barometric law; i.e.: 


n(h) = n(h ) 
o 


T(h ) 

x 2_ 

T(h) 



-3 

cm 


(2.23) 


o 

where n(h) and n(h Q ) are the concentrations at heights h and h^, re- 
spectively; T(h) and T(h Q ) are the corresponding neutral temperatures 
in °K and H(h) is the scale height. 

The vertical distribution of temperature, T(h) , can be expressed 
in terms of the exospheric temperature using the empirical formula 
(Bates[l959], Jacchia[l965]): 


T(h) = T - (T - T ) exp[- a(h - h )] °K (2.24) 

® O O 

where T = T(h ). T is the exospheric temperature, h and h are ex- 
o o 00 - o 

pressed in kilometers and O' is a parameter with dimensions km , differ- 
ent for each profile, the value of which is: 

*2 

a = 0.0291 exp(- 
T ro - 800 

where 5 = 7 2 (with in °K) 

750 + 1.722 X 10~ 4 (T - 800) 2 

CO 

On the other hand, the exospheric temperature can be expressed in terms 
of only the solar input and local time using the formula, (Nisbet[l97l]): 
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T_ = 


£ 

j=o 


COS 


2?rjHL 

24 


)*£ 


sm 


/ 2rtjHL \ 
' 24 ' 


(2.25) 


j=l 


where C . and S . are coefficients given in terms of the daily average of 
*3 J 

the 10.7 cm solar flux, and HL is the local time. 

The validity of the empirical formula (2.24), which gives the neu- 
tral temperature profile, was checked by especializing (2.25) for 
HL = 12 and for three different levels of solar activity (daily averages 
of solar flux considered were S = 75 units, S = 150 units, and S = 250 
units) . With the three exospheric temperatures obtained and T q = 355 °K 
it was possible to draw the corresponding profiles. Fig. (2.3) repre- 
sents the profiles obtained with the empirical formula (dashed lines) 
and those derived from the standard model CIRA [l965] (continuous lines) 
which gives all atmospheric quantities in tabular form. It can be 
seen that the agreement between every pair of curves is quite good in 
the high solar activity case and acceptable in the other two. 

Substituting Equation (2.24) into (2.23) this Equation can then be 
written; 


n(h ) X T(h ) r r 

o o m I 

-T )exp[-a(h-h )] exp | k f 

O O L ■/, 


n(h) = 


T - (T 
00 00 


g(h) dh 

T - (T -T )esp[-Of(h-h )] 


oo oo o 


(2.26) 


where H(h) = 


kT(h) 

mg(h) 


Calling the radius of Earth, then the geopotential height above 
the lower boundary of the F2 region is, by definition; 
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T(K) 


Fig. 2.3. Vertical distribution of temperature at noon for low solar 
activity (S = 75 units), medium solar activity (S =150 units), and 
high solar activity (S = 250 units). Dashed lines; The profile was 
obtained by using the empirical formula given by Bates[l959], Jacchia 
[1965]. Continuous lines: The profile was derived from the standard 

model CIRA [1965] which gives all atmospheric quantities in tabular 
form. 
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o 


Assuming further, that: 


g(h) = g(h Q ) 


the integral term of (2.26) can be easily solved and the neutral con- 
centration will be given in terms of only the solar input, local time, 
and height; i.e.: 


n(h*) = n F P exp(ah ) 
o 


(2.27) 


where F 


T exp (ah*) - (T - T ) 

CO CO Q 

_ 

o 


and p = 1 + 


aH 


kT 

00 

with H = , the exospheric scale height, 

00 mg 

o 


The linear loss coefficient, j3, is given by the formula: 


e> = Y CN 2 ] ♦ Y 0 [0 2 ] sec 

Ct £ 


(2.28) 


where, according to Hintereger et al.[l965], 
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^N 2 = 2 X 10 ^ cm^ - sec ^ 


-11 3 -1 

y = 3 X 10 cm - sec 

°2 

and the N and 0 concentrations, namely [N ] and [O ] in Equation 
(2.28), are obtained from (2.27). 

The electron concentration profile was derived assuming diffusive 
equilibrium in the ionosphere; i.e.: 


V* (Nv) = 0 


Ignoring the action of horizontal pressure gradients, drifts 
caused by electric fields and assuming that the vertical gradients of 
Nv greatly exceed the horizontal gradients, the above equation takes 
the form: ^ 


5h <N V * 0 


(2.29) 


where is the vertical diffusion velocity, given by: 


W. 


n / 1 dN \ 

D " p\ 2H + N dh ' 


The diffusion coefficient, 
be written as: 



increases exponentially upwards and may 
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D e 
o 


D = 
P 


h-h 

o 

H 


Substituting D and W_ into (2.29) and performing the operation indi- 
p D 

cated by the derivative sign, that equation becomes the classical one 
derived by Ferraro in 1945: 



3_ dN _N_ 
2H dh + 2H 2 


0 


The general solution of this differential equation, for the case of H 
independent of height, is: 


N(h) = A exp[ (h - h Q )/H] + B exp[(h - h o )/2H] 


The two boundary conditions imposed to this solution, in order to deter- 
mine the integration constants A and B, are: 

At h = h N(h ) = 0 

o o 

At h = h N(h ) = N 

max max max 

where h is the height of the peak of the F2 layer, 
max 

After expressing the two integration constants in terms of the boundary 
conditions, the electron concentration is given by: 


N(h) = N 

max 


exp[-(h - h )/2H] -exp[-(2h - h 
max ' _ o 

1 - exp[-(h - h V2H] 
max o 7 


h 

max 


)/2H] 


(2.30) 
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The height of the peak, h max > * s determined at any time of the day 
by the criterion given by Rishbeth and Barron [l960] ; i.e. the height of 
the peak under diffusive equilibrium conditions occurs where: 


— = 0.6 (During daytime) 


and where 


= 0.13 (During the night) 


d is the so-called diffusion rate and is given by: 


D sin 2 Y , 
, p -1 

d = — — - — sec 

H 2 


where D^, the plasma diffusion coefficient, is given in terms of the 

_3 

neutral temperature in °K, and composition in cm by the empirical 
formula: 


D = 1.2 X 10 17 T°- 7 /([0] + [ n 2 1 ) cm2 x sec 


Y is the dip angle and H the scale height. Therefore, h , and conse- 

’ max 

quently the electron concentration as given by (2.30) is expressed in 
terms of the solar input and local time. 

Substituting (2,28) and (2.30) into (2.22) this equation can now 
be expressed in terms of only the solar input. A computer program was 
prepared to find the corresponding daily values of (0) for an assumed 
solar input, such as: 
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(2.31) 


S = S + S_ cosil) t 
o 1 s 

The daily average of the integrated production rate was assumed to be 
proportional to the daily average of the solar flux and therefore: 

Q = Q + Q. coso) t 
o 1 s 

By taking in all calculations realistic values of S and S • i.e. of 

o 1 

D , a set of 32 consecutive values of S was obtained from (2.31) for 
Q 

three different levels of solar activitv. 

We note that the expression given by Nisbet for the diurnal varia- 
tion of the exospheric temperature in terms of local time and daily 
average of the solar radio flux does not include any phase information. 
However, it has been established that there is a systematic one day 
delay in the response of neutral temperature to changes in the solar 
radio flux ( Bourdeau, Chandra, and Neupert [1964] , Chandra and 
Krishnamurthy [ 1967 ] , Thomas and Ching [l969], Herman and Chandra [l969]) . 
Consequently, we want to introduce back delay in the Nisbet formula. 

If we consider again the expression of the exospheric temperature: 

j=o j=l 

its daily average will be: 

T = C (S) 

CO o 
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and taking account of the 1 day delay between T ffl and S, this expression 
can be written, after considering (2.31): 


T = C [S(t - 1)] (with t in days) 

co O 

The diurnal variation of the exospheric temperature can then be 
written: 

+ 2 l c j + si »( 2 fr k )l < 2 - 32 > 

j=l 


Every consecutive value of S was substituted into (2.22) from which the 
corresponding daily average of the weighted mean recombination 
coefficient, (3) was determined. The result was a set of 32 values of 


(3) which were passed through a mathematical best-fit filter from which 
and S^g^were obtained. 

With the purpose of finding realistic values for the coefficients 
of Equation (2.31) under low solar activity conditions, a plot was made 
of the solar radio flux at 2.8 GHz for the period of 6 May 1965 through 
1 June 1965. (Fig. 2.4) The choice of this particular period of time 
was dictated by the quasi-sinusoidal character of the solar flux, i.e. 
by the absence of strong erratic changes, which, in the context of the 
present work, are considered as noise. Daily average values of the 
radio flux were plotted together with the corresponding best-fit sinu- 
soidal curve from which the following values were derived; 
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81 units 


and therefore: 


12 units 



0.15 



Fig. 2.4. Graph of the solar radio flux at 2.8 GHz for the period of 
6 May 1965 through 1 June 1965, and its best-fit curve from which 
the long term average, S Q , and amplitude, were derived. 

Fig. (2.5, Curve a) represents the variation of the daily average 
of the integrated production rate. This curve, drawn conventionally 
with zero phase, was derived using the previously mentioned assumption 
of considering the daily average of the integrated production rate to 
be proportional to the daily average of solar flux, that is: 
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Fig. 2.5. Curve (a). Assumed variation of the daily average of the 
integrated production rate. 

Fig. 2.5. Curve (b). Variation of the daily average of (P) throughout 
32 consecutive days and the corresponding best-fit curve. 

Fig. 2.5. Curve (c). Theoretical variation of the daily average of the 
electron content, derived by assuming a Q = 10^-2 cm — 2 ^ throughout 32 
consecutive days. 
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Curve (2.5b) represents the daily average of (3), the weighted 
mean of the recombination coefficient, and its best-fit from which it 
was found that: 


D^= 0.082 


Since a one day delay between S and the atmospheric temperature was 
assumed: 


6 <s>~ 1 day 


12 -2 

Curve (2.5c) was drawn by assuming a Q q = 10 el. cm , which is a 

reasonable value of the integrated production rate^for the level of 

solar activity considered. By substituting the values of D^, and 

§ / \found in Equation (2.21) the delay of the electron content with re- 
\°) 

spect to the solar input was determined. The result is: 


6 = 1.15 days 


In addition, since 6^ has been derived from an equation of the form: 


<P> x i = Q 
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the sensitivity of the electron content, that is — , will be given by: 

I 

AI _ AQ _ A7P7 

i Q W 

°r: D I = D Q " D ( & > < 2 ‘ 33) 

which gives, in this case: 


D t = 0.068 

The same procedure was applied to a set of values of solar flux 
covering the whole range of solar activity, (from = 80 units to 
S Q = 280 units). In order to facilitate comparison, identical values 
of D q were used for all solar activity levels. The results are synthe- 

/N 

sized in Table 1. 

From the observation of this table the following conclusions can 
be drawn: 

a) The delay of the electron content decreases when solar activity in- 
creases. Using the information given in the first and third columns, 
a plot of the delay 6^ versus the solar flux averaged over several solar 
rotations, S Q , was made (Fig. (2.6)). A best-fit exponential curve was 
passed through all 6^ vs. points which allows for writing the follow- 
ing relationship: 

6^ = 3.6 exp( -0.012 X S Q ) days (2.34) 
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TABLE 1 


s o 

D = 0.15 
Q 

D s 

6 0 .= 1 day 

6 j(days) 

D I 

80 

0.085 

1.2 

0.065 

100 

0.077 

0.990 

0.073 

120 

0.0705 

0.836 

0.080 

140 

0.0671 

0.767 

0.083 

160 

0.0541 

0.570 

0.096 

180 

0.0536 

0.550 

0.0964 

200 

0.0502 

0.490 

0.0998 

220 

0.0345 

0.288 

0.115 

240 

0.025 

0.197 

0.125 

260 

0.0189 

0.145 

0.1311 

280 

0.0148 

0.1 

0.135 


The preceding equation, although representing an empirical fit of 

calculated points, is more convenient than ( 2 . 21 ), not only because it 

is simpler, but also because it expresses 6 ^ in terms of the averaged 

solar flux over several solar rotations only. It should be noted that 

Equation (2.34) is applicable also to cases in which different values 

of D are used, 

y 

b) The sensitivity of the electron content is smaller than the sensi- 
tivity of the integrated production rate and increases with solar 

activity. This is caused by the fact that D, > in Equation (2.33) has 

W 

the same sign as does D (i.e. the recombination coefficient 

Q . 

increases when solar radiation increases) and therefore: 


D < 

I Q 
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DELAY OF THE ELECTRON CONTENT 



LONG TERM AVERAGE OF THE SOLAR FLUX 


Fig. 2.6. Scatter diagram of electron content delay versus long term 
average of the solar flux, and the best-fit curve through the scatter- 
ed points. 
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Thus the variations in ^(3^> moderate the changes in I due to varying 
solar inputs. 

The fact that the daily average of the recombination coefficient 
increases with solar radiation can be understood in the following way: 
When solar activity increases the concentration of increases at all 
levels, producing an increase in {3, the linear loss coefficient. How- 
ever, the height of the peak does not vary substantially. As a matter 
of fact, it was found, after applying the Rishbeth criterion to deter- 
mine the height of the peak to CIRA models 2 and 3, that a variation of 
25% in solar flux produces a variation in the height of the peak of 
only 6.4%. Since the weighting factor in the expression of ((3), i.e., 
the electron concentration, increases in the same amount in the 
numerator and in the denominator does not affect the value of this 
parameter because the height of the peak is practically constant and 
therefore ^3 J will increase as a consequency of the increase in 0. 

As shown in Chapter 5 where experimental observations are examined, 
frequently conclusion (b) is violated, i.e., frequently the relative 
variation in electron content is larger than the relative variation in 
solar flux. This brings out one of the limitations of the simple model 
developed in this Chapter. In Chapters 3 and 4 we introduce the effect 
of neutral winds on the ionization and find that the sensitivities of 
electron content can become larger than those of the solar flux. 
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CHAPTER 3 


THE NEUTRAL WIND FEEDBACK MECHANISM 

In the preceding chapter the ionospheric response to the 27-day 
solar ionizing radiation was determined for the hypothetical case in 
which the action of neutral winds was disregarded. 

The general effect of the neutral wind on the ionospheric F2-layer 
(as discussed by King and Kohl [l965]) is to move the ionization along 
the magnetic field lines. During daytime the meridional wind blows 
poleward and tends to push the ionization down, while during the night, 
the meridional wind blows equatorward and tends to lift the ionization. 
The effect of the downward drift, during the day, is to move the F2- 
layer into a region of larger loss rate and thus to decrease the elec- 
tron concentration and the columnar electron content, the reverse taking 
place at night. 

The above effect of the neutral winds is considered in this chapter 

V.- 

in an investigation of the interaction between the 27-day solar ionizing 
radiation and the electron content. An analytical study of the dynamic 
properties of the upper atmosphere is complicated by the need of solving 
a large number of coupled differential equations ( Stubbe [l970]) . In the 
F-region range of heights the neutral air is a mixture of mainly diatomic 
nitrogen and monoatomic oxygen while the ionosphere contains electrons 
and monoatomic oxygen ions. For each of these particle species one has 
to solve a continuity equation, three components of the equation of 
motion, and the heat-transfer equation. In other words, a set of twenty 
coupled differential equations must be dealt with. 
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The large amount of computer time that would be required to solve 
these equations, the mathematical acrobatics associated with the sta- 
bility of the solution, and finally the uncertainties in the determina- 
tion of the boundary conditions that must be imposed on the problem, 
suggest a "non-standard" procedure to solve this particular problem, 
one which will make use of all possible simplifications without compro- 
mising the result. 

Since the output to the procedure must be the daily averages of 
the electron content throughout the period considered, the main simpli- 
fication consists of choosing from among the instantaneous values of 
the diurnal variation of the electron content, one that can be used as 
an index of the daily average. 

With that purpose, a study was made of the cross-correlation 
coefficient between values of the electron content taken every 30 minutes 
and the corresponding daily average. The data used correspond to the 
Stanford ATS-1 station, where the system used for recording electron 
content data takes one sample every five minutes or 288 samples per day. 
The period chosen was the whole year of 1969. It was found that the 
electron content measured at any hour during daytime correlates fairly 
well with the corresponding daily average (correlation coefficient 
p s 0.7). In particular the electron content measured near the middle 
of the day presented a correlation coefficient of the order of 0.9 when 
compared with daily averages. We chose to work with the content observed 
near 1400 L.T. because at that time the continuity equation can be solved 
for the homogeneous case, that is, the ionosphere can be considered 
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under steady-state conditions. 


Let us now consider a steady-state ionosphere at 1400 L.T. sub- 
jected to a certain ionizing flux. To reduce the number of equations 
to be considered, a model for the neutral atmosphere can be adopted. 
Values of the electron concentration N(h) and of neutral .wind velocities 
V n (h) as well as the resulting electron content can be computed in a 
self-consistent way. If now the ionizing radiation is allowed to suffer 
a step increase the following chain of events is triggered: production 

and consequently ionization goes up immediately. This causes an in- 
crease in ion-drag and a reduction in wind velocity. Since the midday 
wind is poleward, as stated above, and forces the ionization into lower 
heights, a reduction in this wind will allow the plasma to rise causing 
the average recombination coefficient to become smaller. This, in turn, 
allows the electron concentration to increase further and reduces the 
wind velocity even more. There is, therefore, under these conditions, 
a positive feedback in the ionosphere which produces a relative increase 
in electron content larger than the relative increase in ionizing radia- 
tion. After a certain delay, generally accepted to be of about 1 day 
as pointed out in the preceding chapter, the augmented energy input to 
the higher atmosphere will result in a neutral temperature elevation 
which causes a decrease in the pressure gradient forces that drive the 
winds, and consequently, a new reduction in the wind velocity. There- 
fore, the electron concentration will increase even further, increasing 
the ion-drag force so that the wind velocity will again decrease. 

When one considers the positive feedback described above, one is 
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confronted with the question of why the system does not run away. An 


idea of the mechanism that limits the increase in electron content can 
be obtained by considering a much simplified model in which diffusion 
is neglected and losses are simply equated to production. Although 
patently unrealistic, such a model allows one to visualize the limita- 
tion mechanism. One can, under such circumstances, write; 


V 


ax 


= 3 


max 


X N 


max 


(3.1) 


where the subscript serves to indicate the fact that the equation is to 
be applied to quantities at the height of maximum ionization. 





Fig. 3.1. Production rate profiles, corresponding to two different 
levels of solar radiation, and S2, used in the text to visualize 
the limiting mechanism of the increase in electron content when the 
solar radiation changes suddenly from to S 2 . 
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Let us now inquire what happens when the solar radiation changes 
abruptly from a value S to a larger value S . Figure (3.1) represents 

1 <o 

two possible production rate profiles corresponding to these two levels 
of radiation. Under solar activity conditions the height of the 
peak occurs at the point A, where according to (3.1): 


(A = (A 
Snax p max 


X N 


(A 


max 


If now the solar activity suddenly changes to the S 2 value, production 
and ionization will change instantly in such a way that the new situa- 
tion will be represented by the point B where: 


(B = 0 <B 

Tnax p max 


X N 


(B 

max 


(3.2) 


in which 0 ^ = 0 ^ . 

max max 

If neutral winds are not considered,' Equation (3.2) will represent 

the new equilibrium situation. If the neutral wind effect is taken 

into account the variation in N will cause a variation in neutral 

max 

winds velocity which in turn will lead to a variation in the height of 

the peak. Therefore, the point B of Figure (3,1) will move up along the 

curve corresponding .to the solar radiation S . 

z 

At any h corresponding to the range of variation of a in 
max inax 

the interval BC, the production rate is such that: 
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Therefore: 


Ah £ H X In 
max 


(B 

Snax 

(A 

Snax 


(3.3) 


Expression (3.3) indicates that the height increase will stop before 

h reaches the value corresponding to the point C, i.e., the new 
max 

equilibrium will occur when the layer has risen to such a height that 
the new q^ is slightly larger than the q max prior to the increase in 
radiation. 

The stability of the system can be demonstrated in a more convinc- 
ing manner through the use of the computer program described in detail 
in Chapter 4. This program models the reaction of the ionosphere to 
sudden changes in ionizing radiation by solving the continuity equation 
under the action of a given neutral wind and then computing the wind 
velocity consistent with the calculated ionization. The new wind infor- 
mation is used to recompute the ionization, from which a correction to 
the wind velocity is obtained. The iterations are carried out until a 
stable, self-consistent solution is obtained. The program, in fact, 
produces a sequence of solutions that carry q^^ along the q-profile 
from B to C', the new equilibrium point. 

Figures (3.2a) through (3.2c) show three typical situations. These 
curves are plots of production rates versus height for low, medium, and 
high solar activity. Point A represents the production rate at the 
height of the peak under solar activity conditions, while point C 
represents the limiting situation predicted by Equation (3.3). In 
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> 11 ) 


-S, =100 UNITS 

-S 2 = 125 UNITS (with the same composition os in 
the S| case) 

-S 2 * 125 UNITS 


I I I I I 


S, =125 UNITS 

-S 2 =175 UNITS (with the same composition 
as in the S, case) 

-S 2 =175 UNITS 


I I I 


-S, =200 UNITS 

S 2 =250 UNITS (with the same composition 
as in the S| case) 

^Sj =250 UNITS 


0 1 23456789 10 

PRODUCTION RATE 


Fig. 3.2. Plots of production rates versus height for low, medium, and 
high solar activity, showing the equilibrium situations. 
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other words, when there is an increase in solar activity from to S 
Equation (3.3) predicts a new value of production rate at the peak of 
the layer lying between points B and C. The new equilibrium situation 
found by the computer program is described by the point C' showing that, 
as expected, the height of the peak increased substantially; however, 
the value of the production rate at this height has changed only a little 
from its original value at A. If now one allows the neutral atmosphere 
to heat up under the influence of the augmented solar radiation, the 
gas density will increase and so will the production rate at any height. 
This is represented by the curve labeled S' . The increase in neutral 
temperature will reduce the pressure-gradient force and, as explained, 
there will be a further decrease of the wind with consequent additional 
rise in the layer: a new equilibrium point is found at C" . 

Rishbeth[op.cit . 1972] showed that, in absence of wind, the ratio 

P/d at the height of the peak is equal to approximately 0.6 during the 

daytime. When the action of neutral winds is taken into account, this 

ratio changes greatly and becomes a function of the solar activity. 

0 

Figure (3.3) shows the values of ^ found at the height of the peak 

d 

for every solar flux considered in Figures (3.2). It can be seen that 
when solar activity increases, that ratio decreases more and more, 
possibly tending towards the value of 0.6 given by Rishbeth. A best- 
fit curve was passed through all points in such a way to make it assymp- 
totic to the value 0.6, i.e., for a very large solar flux, the wind 
effect is sufficiently small so that the value of P/d given by Rishbeth 
is applicable. 
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Fig, 3.3. Scatter diagram of 0/d versus long-term average of the solar 
flux, and the corresponding best-fit curve through the scattered 
points which allows to define the level of occurrence of the height 
of the peak when the effect of neutral winds on the F2 layer is taken 
into consideration. 


The equation of the best-fit curve is: 


| = 36exp (-0.0115 X S) + 0.6 


(3.4) 


This expression is valid only at the middle of the day. 

The block diagram of Figure (3.4) presents in a schematic form, 
the feedback system just described. Box C receives as inputs informa- 
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tion about solar radiation, neutral temperatures (from which a neutral 
atmosphere concentration profile is calculated) and neutral wind pro- 
files. The output of this box is data on ionization in the form of an 
N(h) profile which will be integrated between the two values of height 
that constitute the lower and upper boundaries of the F2-region to 
give the columnar electron content. The N(h) profile will also serve 
as an input to box B together with the neutral temperatures, used here 
to calculate pressure gradient forces; the output is wind velocity which 
completes the feedback loop to box C. Box A modifies the neutral temper- 


atures in response to changes in solar radiation introducing a suitable 
delay in the process. 



Fig. 3.4. Block diagram of the feedback system described in the text. 

The response time of the ionosphere to the solar ionizing radia- 
tion is defined in the context of this work as the lag between the 27- 
day fluctuations in the daily average of electron content and the 27-day 
variations of solar decimeter flux. 

Since the system of Figure (3.4), which represents the feedback 
mechanism described, provides a simple way for finding the output 
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corresponding to step functions, a sinusoidal solar input of 27-days 
period as that given by Equation (2.31), i.e.: 

S = S. + S, cosii) t 
0 1 s 

will be substituted by its quantized form, that is, by a series of step 
functions . 

The smallest quantizing interval that we can take is 1 day since 
this time is approximately the reaction time of the neutral atmosphere 
to the 27-day solar input. Therefore, we have decomposed the sinu- 
soidal solar input into 27 step functions which are shown in Figure (3.5). 



Fig. 3.5. Quantized form of the sinusoidal solar input taken at the 
input of the system represented in Figure (3.4). 

In order to solve numerically the system represented in Figure (3.4) 
for the quantized solar flux applied at the input it will be useful to 
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express all known functions and parameters involved in the equations 
representing boxes A, B, and C, in terms of that input. 

a) Horizontal components of the pressure-gradient force . 

By definition the ^-component of the pressure-gradient force per 
unit mass is: 


F 


? 



(3.5) 


where and p^ are the density and pressure, respectively of the 

neutral atmosphere, and 
§ is any horizontal coordinate. 

Taking the partial derivative of p n from the perfect gas law: 


p = n k T 
n 


where n is the neutral concentration, 
k is the Boltzmann constant, and 
T is the neutral temperature, 
Equation (3.5) will give: 


F 


§ 




n is given by (2.23), i.e.: 


(3.6) 


SEL-73-023 


50 



TOO 

n(h) = n(h Q ) X - exp 


r 


dh 

H(h) 


Differentiating both sides of the above equation with respect to the 

independent variable § assuming that both n(h ) and T(h ) are indepen- 

o o 

dent of § and substituting it into (3.6) that equation is transformed 
into: 


P 5 = 


- T(h) 


f h 

-j 

{ tV, X 

IjTCh) 


dh 


(3.7) 


Since T(h) is given by (2.24), i.e.: 


T(h) = - (T^ - T q ) exp - ar(h - h Q )J 


the derivative term of Equation (3.7) can be substituted by a more ex- 
plicit form and that expression becomes: 


= -T(h) X 


^co 

as J 


f (h ,T ) dh 


(3.8) 


where f(h,T oo ) is given by the formula: 


f<h,T ) = 


CO' 2 

T (h) 


1 - exp[-a(h - h )] + 
o 


750-1.722 X 10~ 4 (T -800) 2 


+ (T - T )(h - h ) 


[750+1.722 X 10 (T -800) ]“ 


rPli 
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Equation (3.8) represents a general expression for any horizontal 
component of the pressure-gradient force per unit mass as a function of 
exospheric temperature, and height. 

In order to especialize Equation (3.8) for the case of determining 
the two horizontal components of the pressure-gradient force correspond- 
ing to an ionospheric point P (see Fig. (3,6)) defined by the distance 
between the point and the center of the earth, its latitude § and its 
longitude X we choose a coordinate system in which the x-positive is 
directed toward the geographic south, the y-positive, toward the geo- 
graphic east and z-positive vertically upwards. In addition, the 
coordinate x is measured along the meridian, the coordinate y is 
measured along the parallel and the coordinate z radially. Therefore: 

x = - (R e + h)$ 

y = X(R,-, + h) cos 5 
E 

z = (R + h) 

E 

where is the radius of the Earth, and 

h the height over the surface of the earth. 

On the other hand, if we assume that the geographic and geomagnetic 
coordinates coincide and that the frame of reference is fixed in space, 
i.e., it is not rotating with the earth, the longitude X corresponding 
to the point P will be given by: 
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Fig. 3.6. Coordinate system used to define a given ionospheric point 
P, in terms of the distance between the point and the center of the 
earth (R E th) , its latitude i, and its longitude X. 
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where is the rotation frequency of the earth, and therefore, the time 
t and the y coordinate will be interchangeable; i.e.: 


d_ _ 3_/St\ 

Sy ~ 3t \ 3y ) 

a i 3_ 

3y ~ (R e + h)^ cos4 3t 


similarly; 


a x i. 

3x R + h 34 

E 


Therefore, the two horizontal components of the pressure-gradient force 
expressed in terms of the exospheric temperature and height are: 


F 

x 


T(h) /dM f 
r e + h \ at / J 


f(h,T oo ) dh 


(3.9) 


F 

y 


T(h) 

(R„ + h)m cos4 
E E 



fth.T^) dh 


(3.10) 


In order to express the two abo^e equations in terms of the solar 
flux only, a model giving the latitude and local-time variations of the 
exospheric temperature was used. Such a model is that of Jacchia[l965] 
which is based on data obtained from observations of satellite motions 
and gives the following formula for the exospheric temperature: 


T m = 

00 CO 


/I . r_. ( -> cos r Tl - sin r 0 r T \ 

(1 + y Sln 0) I 1 + Y 1 — cos — I (3.11) 

min \ 1 + y sin 0 / 
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where T n is the averaged minimum global value of the exospheric 
min 

temperature (function of solar activity); 

0 and 1] are parameters depending on latitude and solar declina- 
tion; 

T is a parameter given in terms of the hour angle of the sun, 
that is, of the local time; and 
y and r are constants, the values of which are 0.3 and 2.5. 

Equation (3.11) allows to express the two horizontal components of 
the pressure-gradient force, as given by (3.9) and (3.10), in terms of 
only the solar flux. 


b) The production model. 


The production of electrons and 0 + ions is caused by the photo- 
ionization action of the solar ultraviolet radiation. The production 
rate at each point of the atmosphere is given by: 


q(h) = [0] X 




(3.12) 


where 


N. 



is the number of bands in which the solar ultraviolet radia- 
tion has been divided; 

is the effective ionization cross-section of the atomic 


oxygen in the band j ; 

$ is the photon flux at level h of the atmosphere, corresponding 
to the band j; and 

[0] is the concentration of monoatomic oxygen. 
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For a practical computation of q(h) , the solar spectrum was divided 
into 16 bands. The photon fluxes and cross-sections for these intervals 
were taken from Garriott, da Rosa, and Davis [l969]. It was assumed 
that because of the rapid recombination of N+ and 0* ions with electrons, 
the ionization of N and 0 moleculars does not contribute to the total 

Ci £ 

production of electrons and ions. However, the absorption of radiation 

by N and (5 moleculars are taken into account. Therefore: 

J 2 2 


$ . = « 
J 


exp 


T.(0) + T (N ) + T . (0 ) 
J J ^ J ^ 


v = the unattenuated photon flux at the top of the atmosphere. 

J 

T.(0) and T.(N ) are the optical depths of 0 and N respectively, along 
J J 2 2 

the path of the solar radiation, and are given by the expressions: 


Tj (O) = 


(0) J n 


(0) ds 


(3.13) 


(a 


Tj( n 2 > = cr (N 2 ) I n(N 2 )ds 


/ 


(3.14) 


' r j (0 2 ) = ^ 


r <v J ■ 


(0 2 )ds 


(3.15) 


( s ( a ( & 

Here, <7. (0) , o. (N ) , amd cr. (O ) are the corresponding absorption 

J J 2 j 2 

cross-sections. 
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c) The loss model. 


The dominant chemical reactions involved in the loss of 0 ions are 


0 + + N -+ N0 + + N with rate coefficient v„ 
2 N, 


0+0 -> 0* + 0 with rate coefficient y 

2 2 0 r 


■N+ J 


The loss rate of 0 ions may be expressed: 


jt(h) = y N [N 2 ] X [0 + ] + y 0 X [0 2 ] X [0 + ] (3.16) 

2 2 


and calling g the linear recombination coefficient, Equation (3.16) 
can be written A(h) = P X [0 + ], where: 


B = V N X [» 2 ] * Y 0 X to 1 
2 2 


d) The diffusion and molecular viscosity coefficients . 

For the diffusion coefficient D we will use the empirical formula 

P 

given in the preceding chapter, i.e.: 

= 1.2 X 10 17 T°- 7 /([0] + [Ng ] j cm 2 - sec 1 

where the neutral concentrations are given by an expression of the form 

(2.27), and therefore, is expressed in terms of only the solar radio 
flux . 


57 


SEL-73-023 



Similarly, the value of p, the coefficient of molecular viscosity, 
is given by: 

p = 4.5 x 10 4 (t/ 1000)^‘ 71 gr. cm 1 sec 1 

(T is given in °K in both expressions) . 
e) Equation of motion for the ionization . 

The momentum transport equation for the 0 + ions is a mathematical 
expression for the sum of all forces acting over a volume element of 0 + 
ions moving with velocity in a gaseous mixture of neutral particles 
and electrons. 

Those forces are due to the pressure gradients of the 0 + ions 
themselves, to collisions between 0 + ions and electrons, and 0 + ions 
and neutral particles, to external electric, magnetic, and gravitational 
fields, and to the Coriolis effect. Therefore, the equation can be 
written: 

- VP + Ne(l + V X B) - mN v . (.V - "v ) - mNu . ( - "v - ) + mNg 

16 6 in n 

+ 2mN(V X UUg) (3.17) 

mass of the 0 + ions, 
concentration of the 0 + ions, 

NkT t 

the electric field 

the geomagnetic flux density (-BcosY, 0, -Bsin¥) with Y the 


dt 


(mNV) = 


where m 
N 


P = 


E 


— ) 

B 
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dip angle, 

n . and u . are the effective collision frequencies between ions 
le in 

and electrons and ions and neutrals, respectively, 

V e = bulk electron velocity, 

V n = bulk neutral particles velocity, 

^ = acceleration due to gravity (0, 0, -g) , 

UU = Earth's rotation frequency, and 
E 

— ) 

V = bulk ion velocity. 

Neglecting the frictional force between 0 + ions and electrons, the 
Coriolis and the acceleration terms because of their small magnitude 
( Rishbeth [op.cit .1972]) Equation (3.17) becomes: 

O = -k V(NTj + Nmg + eNO? + V X B) - Nmu (V - v ) (3.18) 

1 in n 

Equation (3.18) also applies to the electron motion and, neglecting 
small terms which involve the electron mass, one obtains: 

0 = -k VO^T ) - eN(lf + v X B) (3.19) 

e e 

With the purpose of facilitating the solution of the two last 
equations it will be convenient to express them in terms of the 
diffusion velocity, IT , which represents either the ion or electron 
velocity when the neutral air is at rest. This means that under these 
conditions electrons would interact with the ions only via the electric 
field. 
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Therefore, assuming in (3.18) and (3.19) that: 



and: 


v = 0 
n 


we will have: 




k 


mu . 
in 


J V(tfT.) + J V(NT ) 
N l N e 


Defining the plasma temperature, T , the plasma diffusion coefficient, 


D , and its scale height, H^, by the equations; 


T 

P 


T + T. 
e l 

2 


D 

P 


2kT 
P 

mu . 
m 


H 

P 


kT 

P 

mg 


we can express the drift velocity of the plasma, v , by the well-known 
equation: 


= - d p( 


1 

2H 


1_ 

T 


dr 

I 

3Z 


1 9N 
+ N 3Z 


(3.20) 
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Now the solution of the six scalar equations resulting from the two 
vectorial equations (3.18) and (3.19) is carried out in an analytical 
way, giving for the three components of the ion velocity: 


V = v cos 2 '*' + v sin ¥ cos ¥ + V (3.21) 

x nx nz Xx 


V = V, - v cos ¥ (3.22) 

y .Ly nz 


V = v.sin 2 ¥ + v sin ¥ cos ¥ + v sin 2 ¥ + V (3.23) 
z d nx nz Iz 


where V lx’ , and are the three components of the drift velocity 

caused by the electric field. Under quiet conditions, small geomagnetic 
disturbance index, these drifts are small compared with the drift caused 
by the neutral winds as it was pointed out in the preceding chapter 
and therefore, they will be disregarded in further calculations. 

f ) Equation of motion for the neutral air . 

The balance of forces is established between actions due to iner- 
tial plus Coriolis opposing to the summation of the actions due to 
pressure-gradient, ion-drag, gravity, and viscosity. Therefore, the 
equation of motion of the neutral air, after having neglected the grad- 
force due to the collisions between neutral particles and electrons, 
takes the form: 


3v 


n > — > — ^ 


— ^ 

— * LI 2 

g + V 

p n 


Ny . m 

-!=- (v 


v. ) 
n 


(3.24) 
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where p is the density of the neutral atmosphere. 

The first term on the l.h.s. of the above equation will be neglected 
in the steady-state solution used throughout this investigation; besides 
its small magnitude compared with that of the Coriolis term, as dis- 
cussed by Geisler [op.cit . 1966] and Kohl and King [l967], makes it 
negligible even for time-varying solutions of that equation. 

In addition, the vertical component of the vectorial Equation (3.24) 
virtually reduces to the hydrostatic equation F^ = g because the verti- 

_3 

cal component of the Coriolis term is of order 10 X g and the other 
terms are smaller still ( Rishbeth [op.cit . 1972]). Therefore, only the 
two horizontal components of the equation will be considered. 

Calling f^ = 2111^ sin$ , the Coriolis parameter, and = KN, the 
ion-drag parameter, the two scalar equations take the form: 


a 2 v 

i± — ™ f v 

P 3z 2 c ny 


f v sin^Y + F 
d nx x 


= 0 (3.25) 


o v 

j± ny 

P az 2 


fv — f v + F = 0 
c nx d ny y 


(3.26) 


which must be solved for v and v 

nx ny 

g) Equation of continuity. 


The continuity equation, as it was analyzed in the preceding chapter 
can be written, under steady-state conditions, and assuming that the 
spatial variations of the derivative term are much more gradual in the 
horizontal plane than in the vertical: 
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Calling: 


, = BN t |^(NV z ) 

W, = v, sin^Y 
d d 

W = v sinY cosY 
n nx 

the vertical drifts caused by diffusion and neutral winds, respectively 
and: 



which, together with (3.25) and (3.26) will be solved in the next 
chapter for the electron concentration profile, N(h) , and the wind 
velocity v n (h) . 
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CHAPTER 4 


NUMERICAL METHOD AND BOUNDARY CONDITIONS 

In the solution of the system discussed in the preceding chapter, 
the equations involved are those labeled (3.20), (3.25), (3.26), and 
(3.27); i.e.: 


d pl 2H T 3Z N 3Z / 
' P P ' 


(3.20) 


*2 

a v 

ji nx 

P 


az 


_ + f V 

2 c ny 


f v sin 2 Y + F 
d nx x 


= 0 


(3.25) 


a v 

H ny 

> 3Z 2 


fv - f v + F = 0 
c nx d ny y 


(3.26) 


q . BN ♦ 4* * w T » 


(3.27) 


Rishbeth and Barronfop.cit .1960] developed a technique to solve 

simultaneously Equations (3.20) and (3.27), i.e,, to solve the above 

problem disregarding the action of the neutral wind. The same technique 

can be applied to the present situation if the viscosity terms are 

neglected. Equations (3.25) and (3,26) are then reduced to algebraic 

rather than differential equations and, by combining them with (3.27) 

v and v can be eliminated leaving a pair of simultaneous first-order 
nx ny 

differential equations. 
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The effect of neglecting the viscosity terms is not important at 
low levels, say up to about 300 km, because at such heights neutral 
winds are largely controlled by the Coroilis and ion-drag forces and the 
viscosity effect is negligible. At greater heights the opposite occurs, 
i.e., the rapid upward decrease of the density p causes the viscosity 
term to dominate the equation of air motion. 

In order to estimate the magnitude of the viscous term let us con- 
sider a range of heights in which this term greatly exceeds all others in 
the equation of motion, which then can be written; 



(4.1) 


where v represents either v or v 

n nx ny 

Since p varies with height much faster than p, we will assume ‘that 
only the former parameter is height dependent. Therefore: 



(4.2) 


The second term on the r.h.s. of Equation (4.2) was assumed to be zero 
(Equation (4.1)), and since the variation of p with height is of the form; 


p « exp 



Equation (4.2) becomes; 
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or: 



and integrating once this last equation, one finds that: 


dv 

n 

— « exp 
oz 



(4.3) 


If the derivative of v^ varies according to (4.3), v^ will vary follow- 
ing the same law, consequently the terms: 



in Equations (3.25) and (3.26) can be replaced by the approximate values 


and: 




v 

nx 


v 

ny 


Since f , 
c 


the Coriolis parameter, is independent of height, 


com- 


67 


SEL-73-023 



paring at every height p/pH 2 with f , the ion-drag parameter, it is 
possible to determine the height at which the viscosity dominates. From 
that height up to the top of the ionosphere, the neutral wind velocity 
is taken independent of height because in that range the velocity grad- 
ient must tend to zero otherwise the viscosity would become overwhelm- 
ingly large. 

In view of the above, the viscosity effect will be considered in- 
directly as a boundary condition imposed on the neutral wind profile ob- 
tained by integrating the following pair of equations: 

w = _ D / + 1_ + 
d a \ 2Hp T p dz N 3z ) 

q = j3N + N-^— ( w . + w ) + (w + w ) 

oz ' d n ) \ d n/ Oz 


(4.4) 

(4.5) 


where: 


D = D sin 2 ^ 
a p 


and: 


w 

n 


f F + f F 
d x c y 

f 2 + f , 2 sin 2 'F 
c d 


sinYcosY 


The computer program for solving Equations (4.4) and (4.5) includes 

a series of subroutines which allow the calculation, as discussed in the 

preceding chapter, of the parameters D , H , T , q, 0, F in terms of 

a P p 3T X dw 

height for every solar input. The derivative terms ^ - and — are 

oz oz 

th 

approximated by the finite difference equivalents; i.e., at the j 
step of integration, these terms are: 
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We set the lower boundary at the height h Q = 120 km and the upper 
boundary at h = 600 km. At the lower boundary we assume that the elec- 
tron concentration N_ and the neutral wind velocity ~v are identically 

0 n o 

equal to zero. The upper boundary condition for ionization is that the 
diffusive flux of ionization vanishes at the top of the ionosphere, i.e. 


|Nw , ) =0 

V d /z=h 


(4.6) 


J 

The height interval between h and h is divided into n steps such 

u ® 


that: 


h -h 

CO 0 

n = — — — + 1 
Ah 


where Ah is the height-step chosen, taken to be 5 km throughout all 
these calculations. 

The upper boundary condition can be satisfied for any arbitrary 
value of the electron concentration by simply setting 
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However, since the solution proceeds downwards through the iono- 
sphere step by step the value of the electron concentration found at the 
lower boundary will not, in general, satisfy the lower boundary condi- 
tion of, 


(N) 


z=h 


0 


0 


The only way in which both boundary conditions can be met is by a 
series of systematic trials. The possible types of behavior which may 
arise are illustrated in Figure (4.1). 



Fig. 4.1. Types of behavior of the computed electron concentration 
distributions N(h) . 

Solution I leads to exaggerated values of N. Solution II also 
leads to excessive values of N, but first it passes through a local max- 
imum of electron concentration. Solution III leads to negative values 
of N. 
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Solutions I and II stand for "condition A" (excessive values of N) 
and Solution III stands for "condition B" (low values of N) . 


CONDITION 
Was 0 

Changes to A 
Changes to B 

Was A 

Remains in A 
Changes to B 


TABLE 2 


DECISION 


(N) 


z=h 


(N) 


z=h 


-» (N) 
*■> (N) 


z=h 

CO 


z=h 


00 


- AN 
+ AN 


(N) 


z=h 

00 


(N) 


z=h 


00 


(N) 


z=h 


(N) 

z=h 


- AN 

+ AN ( AN-»AN/lO) 


Was B 

Changes to A 
Remains in B 


(N) 


z=h 


(N) 


z=h 


<N) h (N) v, + 

z=h z=h 

00 00 


AN ( AN->AN/lO) 
AN 


The first passage through the computer will be called "Condition 0" . 
Boolean comparisons in the computer recognize the corresponding condi- 
tion. Hence, if a condition was the so-called 0 and went to condition 

A, the electron concentration (N) . was excessive and a new attempt 

z=h 

00 

is made with 


(N) -* (N) . - AN 

z=h z=h 

00 00 
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where AN is the increment of electron concentration introduced in the 
program. 

It should be remembered that all parameters involved in Equations 
(4.4) and (4.5) can be expressed as functions of the solar radio flux, 

S, only. Thus, specifying S specifies a complete electron concentration 
profile and the corresponding neutral wind profile. The two equations 
were solved using the procedure described above, applying a collection 
of values of S taken from the expression; 

S = 100 + 15 coso) t (4.7) 

s 

Here, we have chosen a situation in which the long term average value 
of solar flux is 100 units and, in the nomenclature of Chapter 2, the 
parameter D g is taken as equal to 0.15. The choice of these values 
allows the comparison of the results of this Chapter with those of 
Chapter 2. Thirty-one consecutive values of S were used as the quantized 
waveform at the input of the system described in detail in Chapter 3. 

All parameters were calculated for the autumnal equinox. The result is 
shown in Figure (4.2) which represents the above solar input and the 
corresponding electron content with its best-fit. From this it was 
found that the sensitivity and the delay corresponding to the electron 
content are: 


D = 0.5 
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6^ = 2.6 days 



Fig. 4.2. 
best-fit 


Solar input and the corresponding electron content with its 
calculated from the theoretical model described in the text 


The noise which appears in the computed electron content is probably 
a consequence of a constraint introduced into the program, with the pur- 
pose of saving computer time: the electron content computed at every 
step is stored before the iterative process converges fully. Neverthe- 
less, the total amount of computer time necessary is very considerable, 
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suggestive of further simplifications. 

If the system discussed in Chapter 3 is assumed to be linear with 
only minor linearities left, a simpler technique can be used for the 
special purpose of determining the sensitivity and delay corresponding 
to the electron content. The technique consists of decomposing the 
sinusoidal input in a series of rectangle functions, as shown in Figure 
(4.3); every one of which produces at the output of the system a wave- 
form with its center of gravity shifted with respect to the input. Cal 
culating the response corresponding to a given input rectangle function 
the others can be inferred from linearity considerations. 



Fig. 4.3. Series of rectangle functions taken as an approximation of 
the sinusoidal solar flux of 27-days period. 

Calling A i and 6 i the area and shift, respectively, corresponding 

th 

to the output waveform caused by the i input rectangle, the total 


shift, i.e., the delay between the output and the input is defined by: 

6 


£ Vi 




(4.8) 


£ a 4 

i=0 
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Equation (4.8) implies the use of the superposition theorem, i.e., that 
the system considered is nearly linear. 

When this technique was applied to the solar input given by 
Equation (4.7) the sensitivity and delay found corresponding to the 
electron content were: 


D = 0.474 
6^. = 2.72 days 

which suggests that the assumption made about the linearity of the sys- 
tem does not depart much from reality when the system is considered as 
a simulation of the ionosphere under low geomagnetic activity conditions 
and with a sinusoidal solar input. As a matter of fact, it is shown in 
Chapter 5, where experimental observations are examined, that in periods 
when the geomagnetic disturbance index is small the sinusoidal solar 
input of 27-day period causes a sinusoidal behavior in the corresponding 
electron content, with the same period. 

The computer time needed when these assumptions are made represents 
about 32% of the computer time needed when the quantizing technique is 
applied. 

Therefore, it was decided to apply the approximate technique to 
several values of solar flux, with the purpose of covering the range of 
variation of solar activity. The results are displayed in Table 3. 
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TABLE 3 


s 

o 

D = 0,15 
s 

6 i 

D I 

80 

3.46 

0.481 

100 

2.72 

0.474 

150 

1.80 

0.414 

200 

1.38 

0.372 


Using the values of S Q and 6 shown in Table 3, a best-fit exponen- 
tial curve was calculated which allows for writing the following ex- 
pression: 

= 6 exp( -0.0075 X S Q ) days (4.9) 

The comparison of these results with those found in Chapter 2, leads to 
the following conclusions: 

a) The delay of the electron content decreases when solar activity 
increases, whether or not the neutral wind effect is considered. How- 
ever, under low solar activity conditions, when the neutral wind effect 
is more important and therefore the feedback mechanism described more 
powerful, the variation between the delay calculated when the wind effect 
is considered and the delay calculated when the wind effect is neglected 
is much larger than the same variation calculated under high solar ac- 
tivity conditions. 

b) The sensitivity of the electron content, i.e. the ratio between 
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the amplitude of the a,c. component and the d.c. term, is larger for all 
solar activity levels than the equivalent parameter, i.e. D g , corre- 
sponding to the solar flux. Nevertheless, as a consequence of the feed- 
back mechanism which increases the sensitivity of the electron content, 
the conclusion found in Chapter 2 concerning the increase of sensitivity 
with solar activity now is modified in the sense that the sensitivity of 
the electron content decreases when solar activity increases, because 
the feedback mechanism is more powerful under low solar activity condi- 
tions. 
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CHAPTER 5 


OBSERVATIONS 

In this chapter a statistical analysis of data is carried out 
looking for a quantitative measure of the ionospheric response to the 
the 27-day variations of the solar input. Specifically, we are interest- 
ed in determining the delay of the response of the electron content to 
fluctuations in the solar flux as well as the magnitude of the electron 
content fluctuation for a given change in solar flux. 

Since the 27-day effect is not the only information present in 
the data corresponding to either the solar flux or the electron content, 
it is crucial to isolate that effect from others. With that purpose we 
will consider the 27-day effect present in the data as a periodic signal 
and the rest as a random noise. 

The strong fluctuations observed in the day-to-day variations of 
the daily average electron content will be smoothed out by passing the 
data through a mathematical 7-days running mean filter which does not 
introduce any phase changes and causes very small attenuation to the 
27-days period component. 

Calling 'f' and ^ the estimated values of electron content and of 
solar flux, respectively, i.e., the values derived from the best-fit 
made over the actual data, we can represent them by the equations: 


/ l' = I + A_ cos(u t — ) (5.1) 

o I s 1 




S = S + A„ cos(u) t - f ) 
o S s S 


(5.2) 


where I ,A_ :#■ ,-S A„ and #_ are calculated with the least squares 
o I ’ I ’ o’ S S 


/ 
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criterion. Once the phases $ and 4 are determined, the delay of the 

1 D 

electron content is given by: 


6 = -L- ( $ - § ) days (5.3) 

J. U) 1 b 

s 

In addition, we define; 

= ( S - § ) 2 (5.4) 

as the variance of the solar flux S, and 

crj = ( I - t ) 2 (5.5) 

as the variance of the corresponding electron content. 

Further we assume that everyone of the calculated amplitudes, A 

s 

or A , results from the superposition of an amplitude, namely a or a , 

I O I 

due to the 27-day variation, and an additional term, A or A , result- 

R S R I 

ing from the best-fitting over a random set of points equally spaced. 

The magnitude of this additional term is related to the dispersion around 
the mean of the points over which the best-fit was performed, i.e., to 

a s or V 

In order to find the relationship between the amplitude of the 
additional term and the corresponding standard deviation, let us con- 
sider Np consecutive values, equally spaced, of a random variable y such 
that: 


y = x - 'x 


Ts- 


(5.6) 


SEL-73-023 


80 



where x represents a random variable ( say S or I ) and ^ its estimate 
value . 

The variance of y is, by definition: 



(5.7) 


But y = 0 as a consequence of the least squares criterion applied to cal- 
culate the estimate value of x. Therefore, the second term on the r.h.s. 
of Equation (5.7) vanishes. Substituting (5.6) into (5.7) and consider- 
ing that according to either (5.4) or (5.5) 

( 

CT 2 = ( X - it ) 2 

X 


one finds that: 


CT = CT 

y x 


On the other hand, if a sinusoid 


(5.8) 


R cos( w t - $ ) 

s R 


is fitted through the consecutive values of y, the best-fitting con- 
dition with the least -squares criterion implies that the expression: 

N 

2 k - R cos( v - vf 

i=l 

must be a minimum. 

In the preceding equations, R is the random amplitude of the sinu- 

-V- * 

soidal term, and § is the phase of the sinusoidal term. This phase is 
R 
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assumed to be a random variable uniformly distributed between 0 and 2it, 
and statistically independent of R. 

Taking the derivative of the above expression with respect to R and 
equating it to zero, as required by the existence of a minimum, one 
finds that 


R = 


N 

£ y. cos( « s t - * R ) 
i=l 


N 


i=l 


cos ( 0 ) t - § ) 

s R 


The amplitude of the additional term is defined as the r.m.s. of R,i.e.: 



(5.9) 


Calling: 


N 

_ P 


N 


P 


p = E y k cos0 v = L y k x cos6 k 

k=l K “ k=l 


its square mean will be; 

N N 


-2 £ 

k=l £=1 

N 


y k y * Xcos0 k cose^ 


N N 


2 P J 

= E y k xcose k + E E y k x y 4 cos^xcosf^ 


k=l 


(k=A) 


k=l JL=1 
(k ? JL) 
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The second term on the r.h.s. of the last expression vanishes, and 


2 . 

cos 0, 


1 

2 


Therefore: 



(5.10) 


Substituting this result in (5.9), after taking account (5.8), the 
amplitude of the additional term is given by: 


Therefore: 


and 




(5.11) 


(5.12) 


The error distribution in phase can be found with the aid of the 
Figure(5.1) which represents the triangle formed by the three amplitudes 
involved; i.e., the amplitude of the additional term, A taken with a 
random phase 8 assumed to be uniformly distributed between 0 and 2jt, the 
amplitude caused by the 27-day period component present in the data, a^, 
and the resultant of both, the calculated amplitude, . 
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Fig. 5.1, Sketch of the amplitude measured, A , resulting from the geo- 
metrical addition of the amplitude caused by fhe 27-days effect, a , and 

the amplitude of the random term, A . 

K 

The angle 8 represents the error in phase with which A is taken 
x 

at the moment considered in Figure(5.1). It is given by: 



arcsin 


A R sin0 v 

A / 

x ' 


(5.13) • 


And calling: 


A 


x 



the signal-to-noise ratio. 


the mean value of 0 can be written, 

e 

x 



27 

2^2 


\ arcsin 

J 0 



d8 


days 


Therefore, the dispersion in phase will be: 
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(5.15) 


V 


e 2 + e 2 

e T \ e s e i 


The sensitivity of the electron content for a given variation of 
the solar flux is defined: 


(f ) 


(5.16) 


but, according to (5.14), (5.11), and (5.12) the signal-to-noise ratios 
corresponding to the solar flux and to the electron content are: 


Sc 


and: 


g T = 


1 N 

/ p 1 

f A s ) 

/ 2 * 

{ cr I 

F 

v S ' 

nr 

/ A_ \ 

/ p i 


1 2 ’ 

K / 


(5.17) 


(5.18) 


Therefore, eliminating A and A_ from the last three expressions, 

O I 

Equation (5.16) can be written; 


4f)(f) 


(5.19) 


where d and d are called deviation indices, defined by the expressions 

X s 


d, = 


I I 


d s = S 

o 
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A computer program was prepared, which, using data corresponding 
to the year of 1969 for stations at Stanford ATS-1 and ATS-3, Hawaii, 
Cold Bay, Fort Collins, Rosman, Edmonton, and Urbana, calculates for 
every 32 consecutive values of S and I, and with a sliding interval of 
5 days, the set of values given by the expressions (5.1), (5.2), (5.3), 
(5.4), (5.12), (5.17), (5.18), and (5.19) as well as the cross-correla- 
tion function, for different time-shifts, between the solar flux and the 
electron content data. 

The reason why 1969 was chosen is that during that year the 10.7 
cm solar flux exhibited well developed fluctuations with periods of 
about 27-days and with a long-term component fairly constant, as seen 
in Figure (5.2). 

Defining the total signal-to-noise ratio, g T , as the square root 
of the product of the signal-to-noise ratios corresponding to the solar 
flux and to the electron content, i.e., 

g T = ( g s X gj > 1/2 (5.20) 

it was possible to plot this parameter throughout the year period con- 
sidered. 

When the signal-to-noise ratio as defined above, exceeds some 
arbitrarily chosen level one can expect that the phenomena observed are 
dominated by solar flux fluctuations and the measured delays can be com- 
pared with the theoretical ones computed in Chapter 4. 

The major degradations of signal-to-noise ratio appear to be 
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associated with geomagnetic disturbances as can be seen f rom Figure(5 . 3) 
in which the geomagnetic disturbance index, A , is plotted together with 
g^ corresponding to the eight stations mentioned above. It can be seen 
that there is a clear anticorrelation between these quantities. The dis- 
continuities that appear in some g T curves are caused by gaps in the 
data . 

The effect of variations in the geomagnetic disturbance index on 
the ionospheric response can be seen by comparing the plots corresponding 
to the delay and sensitivity functions derived from equations (5.3) and 
(5.19), respectively, against the plot representing the geomagnetic 
disturbance index. 

If the delay derived from Equation (5.3) is compared with that 
associated with the maximum of the correlation coefficient between the 
solar flux and the electron content, one finds that they coincide with 
one another during periods of low geomagnetic activity index. During 
disturbed periods there is a discrepancy between the two delays but they 
follow the same general trend. 

Figures (5.4a) and (5.4b) represent two sets of delays and sensi- 
tivity functions corresponding to Stanford ATS-1 and Hawaii stations, 
respectively. Neither of these stations present any long gap in the data 
and therefore the functions plotted are continuous throughout the period 
considered. They are representative of the behavior of other stations 
and. lead to the conclusion that the actual value of the time-delay 
between the electron content and the 27-days period fluctuations in solar 
activity can be observed only if the solar flux exhibits a well developed 
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fluctuation with 27-days period and if the magnetic activity is low. 

It can be seen in those figures that the delay hovers between 1.8 
and 2.1 days during the period of large signal-to-noise ratio correspond 
ing to the autumnal equinox. 

The long term average of the solar flux during 1969 was close to 
150 units. The delay found from theoretical considerations for that 
solar activity, Table 3 ( Chapter 4 ), was 1.8 days which is in agree- 
ment with the observed values in Figures (5.4a) and (5.4b). 

Thus, from the above observations, it becomes apparent that under 
appropriate conditions it is possible to predict the ionospheric delay 
with respect to the solar forcing function; however, the prediction of 
the ionospheric sensitivity depends on additional factors not included 
in the model developed in this work, as can be seen by an examination 
of Figure (5.5). This figure represents the electron content behavior 
at three different stations during the favorable period of 22 May 1969 
to 29 September 1969. During this period, as we have pointed out above, 
the solar radio noise varied in an almost sinusoidal fashion and four 
27-day cycles are clearly seen. One would expect that the electron con- 
tent fluctuation at the different stations would be proportional to the 
corresponding solar radio noise fluctuations. This however, is not the 
case. The first and largest cycle of solar radio noise produced a rela- 
tively small effect in the ionosphere, especially at Urbana and Hawaii. 
An examination of the fluctuations of the 304 % He II line during this 
period, helps to partially explain the discrepancy observed: It appears 
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Fig. 5.4a and 5.4b. Plots of delays and sensitivities corresponding to Stanford ATS-1 and Hawaii 
stations, during the year of 1969. 








that the radio flux did during the first cycle depart considerably from 
the behavior of the EUV, showing an exagerated amplitude. 

The difference of behavior of the ionosphere at various places, in 
response to a given fluctuation in solar radiation can obviously only be 
explained by differences in local conditions, differences which have not 
been taken into account in the theoretical study made in this work. It 
was shown in Chapter 2 that in the absence of the wind feed-back mecha- 
nism the electron content sensitivity is expected to be less than 1, and 
in the more complete investigation in Chapters 3 and 4, it was shown that 
the inclusion of the feed-back mechanism will cause the sensitivity to 
become larger than unity. The observed sensitivities at Stanford and 
Hawaii during the period in question, are all larger than one suggesting 
that indeed, the wind feed-back mechanism is present. 
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CHAPTER 6 


CONCLUSIONS 

The general purpose of the present work was to examine the dynamic 
interaction between electron content and solar flux, and in particular 
to determine the ionospheric response to the 27-day variations of the 
solar flux, i.e., the sensitivity and delay corresponding to the electron 
content . 

The dependence of the level of solar activity to changes in the 
ionospheric response has been theoretically analyzed for both the case 
in which the action of neutral winds on the F2 layer is disregarded and 
for that in which the wind effect is considered. The comparison of the 
results leads to the conclusion that neutral winds play a dominant role 
in the mentioned dependency. The reason lies in the positive feed-back 
mechanism discussed in Chapter 3. 

A technique has been developed in Chapter 2 for solving the con- 
tinuity equation for the daily average of the electron content when the 
action of neutral winds is disregarded. The technique consists in taking 
daily averages of the continuity equation, written in terms of the elec- 
tron content. This technique leads to two important simplifications: the 
integrated transport term cancels out and the magnitude of the time- 
derivative of the daily average of the electron content is negligible 
when compared with the magnitude of the other terms. From the simplified 
version of the continuity equation it was possible to derive the sensi- 
tivity and delay corresponding to the electron content in terms of the 
solar input and of the daily average of the weighted mean of the re- 
combination coefficient. 
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By adopting a model for the neutral atmosphere and using an em- 
pirical expression relating the neutral temperature profile to the solar 
flux, it was possible to express both delay and sensitivity in terms of 
only the solar input. 

When the action of neutral winds was considered the continuity 
equation of ionization and the equation of motion of the neutral atmo- 
sphere were solved for the ionization and the neutral wind velocity 
profiles in a self-consistent fashion. The simplifications made in both 
equations include: 

a) The electron content measured at 1400 L.T. was taken as an in- 
dex of the daily average of the electron content. This simpli- 
fication allowed considering the continuity equation under 
steady-state conditions. 

b) The viscosity effect on the neutral wind velocity profile was 
considered only as a boundary condition. This simplification 
made it possible to convert the two differential equations of 
the horizontal components of the neutral wind velocity into a 
pair of algebraic equations, the solution of which was carried 
out analytically. 

The technique for integrating the reduced pair of equations result- 
ing from the above simplifications was discussed in detail in Chapter 4. 

With the same model of the neutral atmosphere adopted in Chapter 2 
and the same empirical expression relating the neutral temperature pro- 
file with the solar flux, it was possible to develop an elaborated set 
of formulae which allowed the calculation of all parameters involved in 
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the equations to be integrated in terms of only the solar input. There- 
fore, the only datum required to feed the computer program, prepared 
for integrating the pair of differencial equations, is the solar flux. 

A sinusoidal solar input with 27-days period applied at the in- 
put of the non-linear system proposed in Chapter 3 was quantized, i.e., 
it was considered as a series of step functions. At the output of the 
system one finds a wave-form, representing the electron content , from 
which it was possible to derive the corresponding sensitivity and delay. 
Furthermore, with the purpose of saving computer time, a major simpli- 
fication was done consisting in linearizing the system. However, the 
results did not show much sensitivity to this simplification. 

The delays in the ionospheric response decrease with increasing 
solar activity in' both the "no wind" model and the one in which wind 
effects are considered. In this latter case the delays are somewhat 
larger. 

The "no wind" model predicts a sensitivity smaller than unity 
while the effect of the wind causes this sensitivity to be larger than 
one, i.e., the relative fluctuations in electron content become larger 
than those in solar radio flux. With no winds, an increase in solar ac- 
tivity levels leads to an increase in sensitivity, the opposite occurr- 
ing in the model that takes winds into account. 

The experimental verification of the theoretical results described 
above is difficult in view of the presence of a number of other effects 
not included in the model. Such effects were treated as noise while the 
response of the ionosphere to the fluctuation of the solar flux was 
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taken as the desired signal. The observed data were examined only when 
the " signal-to-noise ratio" was sufficiently high. 

This requires the presence of two conditions: 

a) Solar flux fluctuations with sufficient amplitude. 

b) Low amplitude of the perturbations. 

Condition a) occurs, in general, during geomagnetically quiet periods. 

When the above conditions coexist, the ionosphere responds in 
approximately the manner predicted by the model described in Chapter 4; 

- The observed delays are of the magnitude predicted by the theory. 

- The sensitivities are larger than one. 
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